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Abstract 

We report a detailed study of the electronic structure of the MgB2 with 
different lattice constants by using the full-potential linearized augmented 
plane wave(FPLAPW) method. It is found that the lattice parameters have 
great effect on the a band of Boron. Our results indicate that increasing the 
lattice constant along the c axis will increase the density of states(DOS) at 
the Fermi level, shift the a band upward, and increase the hole number in 
the a band. So, the superconducting transition temperature T c will be raised 
correspondingly. Changing the lattice constant along a axis has the opposite 
effect to that of the c axis. Our result is in agreement with experiment. A 
possible way of searching for higher T c superconductor has been indicated, 
i.e., making MgB2 to have longer c axis and shorter a, b axis by doping. 



PACS Number : 74.25. Jb, 74.70.Ad 
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The discovery of superconductivity in magnesium diboride (MgB 2 ) [1,2] has stimulated 
world wide excitement. The superconducting transition temperature for MgB 2 (T c ~ 39 
K) exceeds by almost two times the record values of Tc for conventional Bl- and A15-type 
intermetallic superconductors (SC) [3], which is by far the highest if we exclude copper 
oxides and C 60 based materials. As distinct from the high-temperature superconductor, 
MgB 2 has an exclusively simple composition and crystal structure [4]. The B isotope shift 
of T c reported by Bud'ko [5] and most other early experimental data [6] suggests conventional 
BCS strong-coupling s-wave electron-phonon (EP) pairing. 

Band structure calculations show that the compound is not only quite ionic with a rea- 
sonable DOS, but also has strong covalent B-B bonding (the bonding-antibonding splitting 
due to in-plane B-B hopping is about 6 eV) and thus exhibits strong electron-phonon in- 
teractions^, 8, 9]. Mg s-states are pushed up by the B p z orbitals and fully donate their 
electrons to the boron-derived conduction bands. Hirsch used a hole-superconductivity 
model to explain the high temperature superconductivity in MgB 2 as driven by undress- 
ing of hole carriers in the planar boron p X J/ orbitals in the negatively charged B~ planes 
[10]. The model by Hirsch predicts a positive pressure effect on T c . However, high pressure 
studies show negative pressure coefficient of T c [ll, 12]. Based on an estimate of phonon 
frequencies and band structure calculations, Kortus et al[8] explain the superconductivity 
in MgB 2 as a result of strong electron-phonon coupling and An and Pickett [7] attributed 
it to the behavior of p^-band holes in negatively charged boron planes. Many authors[7-9] 
emphasize the significant role of metallic B states in the appearance of superconductivity. 
According to the McMillan formula for Tc [13], the high transition temperature is probably 
due to a high density of states at the Fermi level, N(E F ), high averaged electron-ion matrix 
elements, as well as high phonon frequencies, which increase for light elements and depend 
on B — B and M — B bonding. 

Past experience said that T c should be strongly affected by doping or changing the lattice 
constants. Some experiments concerning substitutions on the Mg site have been done [14, 15]. 
Substitution of Al leads, apart from doping with electrons, to a compression of the structure 
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due to the difference in ionic radius between Al and Mg. The compression is anisotropic 
and an Al content of x=0.1 leads to a structural instability. The rate of compression along 
the c axis is about twice as much as along the a axis[14]. Both hole doping by substitution 
of Li + [15], and electron doping by substitution of Al 3+ [14] for Mg 2+ led to a decrease of 
T c . In fact, Al doping will destroy bulk superconductivity when the Al content x is larger 
than 0.3[14]. There are also many experiments concerning Boron sublattice doping [16]. 
The lattice parameter a decreases almost linearly with increasing carbon content x, while 
the c parameter remains unchanged, indicating that carbon is exclusively substituted in 
the Boron honeycomb layer without affecting the interlayer interactions. T c also decreases 
about linearly as a function of the carbon concentration. A recent experiment had also 
found BiB 2 not to be a superconductor [17]. BiB 2 and MgB 2 are isostructural, and their 
valence electron number is the same. So, the basic change is due to structural factors, i.e. 
lattice parameter, a, and interatomic distance, c/a. Recently, high pressure studies show 
that T c will decrease with increasing pressure [11, 12]. Compression will decrease both lattice 
constant a and c, but it is not yet know which one of them is responsible for decreasing T c . 
So, it is very interesting to investigate the Boron band structure, especially the n and a 
band, as a function of lattice constants a and c. 

We have used the highly accurate all-electron full-potential linear augmented plane 
wave method[18, 19]. The standard local density approximation to the electron exchange- 
correlation potential was used together with the generalized gradient corrections of 
Perdew[20]. The muffin-tin sphere radii (R) of 2.00 a.u. and 1.50 a.u. were chosen for 
the Mg and B atoms, respectively, with a cutoff RK ma;E =8.0. 

MgB 2 is isomorphous with A1B 2 [4]; the lattice constants of the hexagonal unit cell are 
ao=3.0834 and Co=3.5213 A. From the existing experiments we can see that high hydrostatic 
pressure has an anisotropic influence along a and c. The compression along the c axis is 64% 
stronger than along the a axis, which is in line with the weaker Mg-B bond[21]. Thermal 
compression along the c axis is about twice the one along the a axis. So, we fix a=3.0834A, 
and vary lattice constant c. The total energy- vs-c relation is shown in Fig.l. It can be seen 
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that c=co=3.5213 A with the lowest total energy agrees with the experimental result. 

The band structure of MgB 2 is shown in Fig. 2. For convenience, we set E/=0. From it we 
can see that there are two kinds of bands: the a and n band. Both of them are contributed 
by Boron. The a band along T-A is double degenerate. Change of lattice constant will not 
change the symmetry, and so will not change its character. We found when c—IAcq, the a 
band has a very small dispersion along T-A. But it will increase with decreasing c. Thus 
increasing c will strengthen the 2D character of the a band, and if c approaches oo, MgB 2 
becomes an ideal 2D B 2 layer. Another evident change is that the a band will shift upward 
with respect to the Ef when we increase c. For c=0.8c , the a band is below the Fermi 
energy, and the a bonding state would be completely filled. For c=0.9c , the a band will 
cross the Fermi energy at the Y point. With increasing c, the a band will shift upward, and 
so will have more holes. On the other hand, increasing c will decrease the hopping coefficient 
t ppa in Kortus's TB model[8] for the ir band, and will cause the ir band dispersion along T-A 
decrease. Comparing with Ef, increasing c will decrease the position of bonding n band at 
M and T point, and will raise the position of bonding n band at A and L point. 

We have also studied the effect of lattice constant a on the electronic structure, and found 
it to be opposite to c. With decreasing a, both bonding and antibonding a band will move 
upward with respect to E/, and the splitting between bonding and antibonding a band will 
increase. When c=0.8c and a=a , the bonding a band is below the Fermi energy. Fixing 
c and decreasing the a, the bonding a band will move upward. From Fig.2(d), we can see 
that, when c=0.8c and a= 0.9a 0) its a band lies higher than it in Fig.2(a), and there are 
holes in the bonding a band. The increase of splitting between bonding and antibonding 
a band is due to increase of the hopping coefficient with decreasing a. The possible reason 
why a band moves upward is electrostatic effect. Being different from a band, the bonding 
7r band will move downward with respect to the Ef when a decreases. Compression will 
decrease both the a and c. But a bond of B-B is much stronger than the Mg-B bond, so, 
it is natural to expect that compression will mainly cause the decrease of c[21]. Thus, we 
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can conclude that compression will shift the a band downward and decrease the number of 
holes in the a band. 

We have also studied the relationship between the density of states at the Fermi level(N/) 
and the lattice constants, which is shown in Fig. 3. It is found that N/ will increase with 
increasing of the lattice constant c, which is in agreement with other theoretical results[22, 
23]. Many theories and experiments show that the MgB 2 is a BCS-like superconductor. 
So the electron-phonon coupling constant A, which enters the BCS equation, is very impor- 
tant [13], and is proportional to N/. Since N/ will decrease for smaller c, T c will be decreased, 
which is in agreement with experiments [11, 12]. The states at the Fermi level, responsible 
for superconductivity, show two different orbital characters: p — a bonding (column-like FS 
around T — A) and p z , which has 7r-bonding and antibonding characters on the basal and 
on the top {k z = ir/c) planes, respectively. Which one of them is the most important for 
the superconductivity is not yet known. In order to answer this question, we calculate the 
partial DOS of n and a band of MgB 2 with different lattice constants c. We find that with 
a larger c, N/(cx) will increase, whereas N/(7r) will decrease. Here, the N/(a) (N/(-7r)) means 
the density of a (it) state at Fermi energy. So, we can say that the a band of Boron plays 
a more important role for superconductivity than the n band. The p CT and p w bands move 
with respect to each other with increasing c, thus inducing charge transfer between these 
two bands. It is the reason why the N/(a) has an opposite trend to N/(7r). Our result is in 
agreement with Goncharov and Bud'ko et al.,'s result [24]. 

In summary, we calculated the electronic structure of MgB 2 for different lattice constants. 
Our results show that increasing lattice constant c will increase the DOS at the Fermi level, 
shift the a band upward, and thus increase the hole number in the a band. Decreasing c 
will lower the T c . The effect of changing a is opposite to c. For constant c, a shorter a will 
raise T c . The main influence of compression is to decrease c, and thus will decrease the T c , 
which is in agreement with the experiments [11, 12]. A search for a higher T c superconductor 
by doping MgB 2 , would require a longer lattice constant c and shorter a, b. A possible way 
would be a suitable substrate, which forces MgB 2 film to have longer c but shorter a, b. 
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Figure Captions 

Figure 1 The total energy of MgB 2 as a function of lattice parameter c/c . c = 3.52 13 A. 

Figure 2 The band structure of MgB 2 with different lattice constants. The circles 
represent the Boron a band. ao= 3.0834A, co= 3.5213A (a) a=%, c=0.8co; (b) a=ao, 
c=co; (c) a=a , c—lAco] (d) a=0.9ao, c=0.8co. 

Figure 3 Density of states near E/ for MgB 2 with different lattice constants, c/c =0.8, 
1.0 and 1.4. Fermi level is set at the zero. 
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